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Abstract 

Background: Acquisition of exogenous genetic material is a key event in bacterial speciation. It seems reasonable to 
assume that recombination of the incoming DNA into genome would be more efficient with higher levels of 
relatedness between the DNA donor and recipient. If so, bacterial speciation would be a smooth process, leading to a 
continuous spectrum of genomic divergence of bacteria, which, however, is not the case as shown by recent findings. 
The goal of this study was todetermine if DNA transfer efficiency is correlated with the levels of sequence identity. 

Results: To compare the relative efficiency of exogenous DNA acquisition among closely related bacteria, we carried out 
phage-mediated transduction and plasmid-mediated transformation in representative Salmonella strains with different 
levels of relatedness. We found that the efficiency was remarkably variable even among genetically almost identical 
bacteria. Although there was a general tendency that more closely related DNA donor-recipient pairs had higher 
transduction efficiency, transformation efficiency exhibited over a thousand times difference among the closely related 
Salmonella strains. 

Conclusion: DNA acquisition efficiency is greatly variable among bacteria that have as high as over 99% identical 
genetic background, suggesting that bacterial speciation involves highly complex processes affected not only by 
whether beneficial exogenous DNA may exist in the environment but also the "readiness" of the bacteria to accept it. 
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Background 

Speciation of bacteria, i.e., their divergence from the 
ancestor and evolution into new species, is facilitated by 
incorporation of laterally acquired genetic material into 
the genome, as demonstrated by comparative studies of 
model bacteria, such as Salmonella and other enteric 
bacteria [1-3]. Bacteria of the Salmonella genus have been 
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frequently used as ideal research models of bacterial 
genomic divergence and evolution owing to several 
advantages that fit such studies, including the close gen- 
etic relatedness and, conversely, distinct pathogenic prop- 
erties of these bacteria. In addition, the extraordinarily 
large number of known genetic and pathogenic types of 
Salmonella makes the comparative studies informative 
and feasible. 

Salmonella diverged from Escherichia coli more than 
100 million years ago [4-6]; to date, more than 2500 
serologically defined types, called serotypes or serovars, 
have been documented [7,8]. Genomic comparisons reveal 
that different Salmonella serotypes have highly similar 
overall genome structures, which are also similar to that 
of E. coli [9-12], validating the long speculated high level 
genomic conservation during long evolutionary times. 
When more detailed comparisons are made systematically 
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over the whole genome, specific differences can be found 
among even very closely related bacterial lineages 
[11,13-16]. Genomic differences can be classified essen- 
tially into two main categories: DNA sequence variations 
and distinct sets of insertions between the bacteria. It is 
well documented that monophyletic Salmonella serotypes 
{versus polyphyletic serotypes that contain genetically 
distinct lineages but are classified together due to their 
common serological features, such as Salmonella paraty- 
phi B containing both single-host and broad-host infecting 
lineages that cause paratyphoid and gastroenteritis, respect- 
ively, in humans) may have "hallmark" insertions, such as 
the 134 kb SPI7 in most S. typhi isolates (S. paratyphi C 
and S. dublin have different versions of SPI7) [17-20]. 
Within a monophyletic Salmonella serotype, individual 
sub-lineages may also have their own unique insertions or 
combinations of them, as exemplified by the DT104 phage 
type of S. typhimurium for their possession of three inser- 
tions ST104, ST104B and ST64B, a combination not seen 
in other S. typhimurium sub-lineages [16,21]. All these facts 
suggest that, although bacteria have many chances to 
contact exogenous DNA in the environment and may 
even possibly have it moved into the cell via mecha- 
nisms such as bacteriophage-mediated transduction and 
plasmid-mediated transformation, the recipient may or 
may not have it incorporated into the genome regardless 
of whether or not the recipient might gain better fitness 
to the environment by having it. For example, although 
in vitro recombination efficiency depends linearly on 
the levels of the sequence similarity of the DNA strands, 
whether this linearity also exists in the exogenous DNA 
incorporation process inside a bacterial cell is unknown. 

In this study, we assessed the efficiency of exogenous 
DNA acquisition by phage-mediated transduction and 
plasmid-mediated transformation in selected bacterial 
strains to establish whether DNA transfer efficiency 
might be correlated with the levels of sequence identity. 
We found that the efficiency was remarkably different 
among even very closely related bacteria. We conclude 
that DNA acquisition efficiency is greatly variable 
among bacteria that have as high as over 99% identical 
genetic background, which implies that bacterial speci- 
ation involves highly complex processes affected not 
only by whether beneficial exogenous DNA may exist in 
the environment but also the "readiness" of the bacteria 
to accept it. 

Results 

Transduction frequency and its correlation with sequence 
similarity: general tendency and strain-specific efficiency 

We transferred five Tn70-inserted genes, including flgL, 
treA, ompC, hisA and purG, from S. typhimurium LT2 
to the recipient strains (Table 1). To assess whether the 
five genes might be representative of the whole genome 



regarding the divergence among the Salmonella strains 
used in this study, we concatenated the sequences of the 
five genes and constructed a phylogenetic tree (Figure 1A) 
for comparison with the tree that was based on the core 
genes of the whole genome (Figure IB). For the strains for 
which no whole genome sequences were available for tree 
construction, we used the published sequences of strains 
of the same serovar, such as S. agona strain SL483 instead 
of SARB1 (Table 1). The two trees demonstrated essen- 
tially the same genetic relatedness among the bacteria, 
justifying the use of these Salmonella strains, although 
the five genes had different levels of divergence among 
the Salmonella strains (Table 2). 

Recombination frequencies as reflected by the relative 
numbers of transductants on LB plates varied among 
thefive genes in a given recipient; however, the differ- 
ences between strains were much greater (Table 3). As 
would be expected, S. typhimurium 14028S exhibited 
the highest frequency of recombination with the DNA 
from S. typhimurium LT2 among all recipient strainsand 
we did not get any transductants from S. bongori SARC12 
(the most distantly related species with S. typhimurium 
among all Salmonella lineages), althoughthis latter strain 
did express the 012 antigen after transformation by 
pPR1347 (See Methods). Notably, however, although 
14028S has nearly identical genomic sequence to LT2, 
its recombination frequency was as low as merely one 
tenth that of LT2. A similar situation was seen with S. 
enteritidis LK5, whose sequence divergence is only 
0.012 from S. typhimurium LT2 but its recombination 
frequency with LT2 was even lower than S. arizonae 
SARC5 with LT2, although S. arizonaeis a serovar from 
subgroup Ilia, which is much more divergent than S. 
enteritidis from S. typhimurium (see Figure 1). Another 
notable observation was that we did not obtain any 
transductants from S. typhi Ty2 even though we tried 
numerous times, which however was consistent with 
our previous findings that we had to disrupt the mutL 
gene of S. typhi Ty2 in order to obtain transductants in 
the transduction experiments with the donor DNA from 
S. typhimurium LT2 when we were constructing a phys- 
ical map for S. typhi Ty2 [18]. 

As a whole in the experiments described above, al- 
though we saw a general tendency of higher transduc- 
tion efficiency with closer genetic relationships between 
the donor and recipient bacterial strains, a linear model 
was not supported by our data (i.e., the relationship was 
not linear; Figure 2). Especially, the results out of S. 
typhimurium LT2 and 14028S seem to demonstrate that 
the transduction efficiency is rather strain-specific. We 
thus considered evaluating bacterial strain specificity by 
involving multiple strains of a Salmonella lineage and 
between very closely related lineages, i.e., S. gallinarum 
and S. pullorum. 
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Table 1 Bacterial strains used in this study* 


Strain 


Accession number 


Relationship with S. typhimurium LT2 


Reference 


S. typhi murium LT2 


AE006468.1 


- 


[38] 


S. typhimurium 14028S 


CP001 363.1 


Same serovar (1,4,[5],12:i:1,2) 




S. typhi Ty2 


AE014613.1 


Same subgroup (I) 


[18] 


S. paratyphi A ATCC91 50 


CP000026.1 


Same subgroup (I) 




S. paratyphi B SPB7 


CP000886.1 


Same serogroup (04) 




S. paratyphi C RKS4594 


CP000857.1 


Same subgroup (I) 




S. agona SL483 


CP001 138.1 


Same serogroup (04) 




S. agona SARB1 


N/A 


Same serogroup (04) 




S. duo//nCT_02021853 


CPO0 1144.1 


Same subgroup (I) 




S. cfub//n 1TT3627 


N/A 


Same subgroup (I) 




S. enteritidis P 125 109 


AM933172.1 


Same subgroup (I) 




S. enteritids LK5 


N/A 


Same subgroup (I) 




S. pullorum RKS5078 


CP003047.1 


Same subgroup (I) 


[23,24] 


S. gallinarum2S7/9] 


AM933173.1 


Same subgroup (I) 


[25] 


S. arizonae RKS2980 


CP000880.1 


Closely related subgroup (Ilia) 




S. arizonae SARC5 


N/A 


Closely related subgroup (Ilia) 




S. bongo/7 NCTC 12419 


FR877557.1 


Distantly related subgroup (V) 




S. bongo/7 SARC 12 


CP006692 


Distantly related subgroup (V) 





*See more detailed information on these bacterial strains at http://www.ucalgary.ca/~kesander. 



Transduction frequency between S. gallinarum and 

5. pullorum: variable efficiency among individual strains 

S. gallinarum and S. pullorum have a common antigenic 
formula, 1,9,12:-:-, the former causing typhoid and the 
latter causing pullorum disease in fowl They were ori- 
ginally treated as separate species [22] but have since 
the mid 1980s been classified into the same serovar of 
the same subspecies (i.e., S. enterica subspecies enterica 
Serovar Gallinarum as separate biovars Gallinarum and 
Pullorum, respectively [8]). However, their distinct bio- 
logical properties, i.e., causing entirely different diseases, 
unambiguously distinguish them as different organisms. 
Recently our work reveals that the two pathogens have ac- 
cumulated distinct sets of mutations, including different 
pseudogenes [15,23], which further demonstrated their 
genetic divergence. We moved the Tn7(Mnserted ompD 
gene first from S. typhimurium LT2 to four wild type 
strains each of S. pullorum (including strain RKS5078 
[23,24]) and S. gallinarum (including strain 287/91 [25],) 
and then from one of the eight strains to the other seven 
strains; this process was repeated for each of the eight 
strains. We observed a tendency that transduction fre- 
quency was higher in a donor-recipient pair of strains 
within S. pullorum or within S. gallinarum than between 
S. pullorum and S. gallinarum strains, although not all of 
these differences were statistically significant (Figure 3A). 
To validate this observation and rule out the possibility 
that a particular genomic DNA segment (e.g., ompD in 



Figure 3A) or some bacterial strains might have biases, we 
used TniO-inserted leu, bio, oxrA and cysA as donor 
DNA, with ompD also included as donor DNA for a 
comparison, in the second set of transduction experi- 
ments with larger numbers of S. pullorum and S. galli- 
narum strains. Again, the transduction frequency was 
lower between S. pullorum and S. gallinarum than within 
S. gallinarum or S. pullorum, but, similar to the above 
experiments, not all of the differences were statistically 
significant (Figure 3B & C). 

Transformation efficiency among individual 
Salmonella strains 

As transduction frequency reflects a combined result of 
DNA taking up capability and recombination efficiency, 
we wondered whether the non-linearity of transduction 
frequency with the level of relatedness of the bacteria 
might have mainly resulted from differential capability of 
the bacteria to take up DNA from the environment into 
the cell. To look into this, we transformed representative 
strains with the plasmid pQE30, which does not have 
homologous sequences with the genomes of the bacterial 
strains used in the study, and compared the transform- 
ation efficiency among them. Intriguingly, the transform- 
ation rates differed more than one thousand times among 
the tested bacteria, e.g., as low as ca. 8.3 x 10" 7 in S. enteri- 
tidis LK5 andas high as ca.1.3 x 10" 3 inS. arizonae SARC5 
(Table 4). 
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Discussion 

Two classes of genomic changes are associated with the 
phylogenetic divergence and evolution of bacteria - 
acquisition of laterally transferred DNA and nucleotide 
substitution, with the former being the primary driver of 
bacterial speciation, i.e., the process and consequence of 
bacterial development into a nascent biological species. 

The biological species is defined on the basis of sexual 
reproduction potential of the organisms [26,27], i.e., each 
species is separated from others through reproductive 
barriers. This definition of species works well for many 
plants and animals but not so much for bacteria owing to 
the mostly non-sexual reproduction modes of bacteria. As 



bacteria acquire much of the genetic novelty through 
homologous recombination of laterally transferred DNA, 
it seems reasonable to presume that recombination effi- 
ciency among bacteria should be falling gradually with in- 
creasing sequence divergence. Several lines of evidence 
exist to support this presumption in Escherichia, Bacillus 
and Streptococcus [28-30] . If so, bacterial speciation would 
be a gradual process, leading to continuous divergence of 
the bacterial genomes without genetic boundaries to de- 
lineate bacteria into discrete clusters, i.e., bacterial species 
would have to be arbitrary. However, our recent work has 
demonstrated the existence of clear-cut genetic boundar- 
ies among bacteria as closely related as the Salmonella 
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Table 2 Sequence divergence of the five genes between 5. typhimurium LT2 and the recipient strains 



Recipient strain 


flgL 


treA 


hisA 


ompC 


purG 


Average 


SD 


S. typhimurium LT2 
















S. typhimurium 14028 s 


0 


0.0006 


0 


0 


0 


0.0001 


0.0003 


S. paratyphi B SPB7 


0.0094 


0.0093 


0.0136 


0.0088 


0.0103 


0.0103 


0.0019 


S. agona SARB1 


0.0105 


0.0146 


0.0163 


0.0202 


0.0085 


0.014 


0.0047 


S. dublin 1TT3627 


0.0084 


0.0088 


0.0081 


0.0088 


0.0082 


0.0085 


0.0003 


S. enteritidis LK5 


0.0126 


0.0088 


0.0122 


0.0202 


0.0082 


0.0124 


0.0048 


S. paratyphi C RKS4594 


0.0063 


0.0128 


0.0108 


0.0202 


0.0113 


0.0123 


0.0051 


S. paratyphi A ATCC91 50 


0.0073 


0.0123 


0.0136 


0.0053 


0.0116 


0.01 


0.0035 


S. fyp/i/ Ty2 


0.0073 


0.0169 


0.0136 


0.0053 


0.01 1 1 


0.0108 


0.0047 


S. arizonae SARC5 


0.0692 


0.0736 


0.0637 


0.0325 


0.0598 


0.0598 


0.0161 


S. bongo/7 SARC 12 


0.1143 


0.0362 


0.061 


0.0554 


0.0769 


0.0688 


0.0293 



lineages [31], indicating that a fundamental question still 
remains regarding molecular mechanisms involved in bac- 
terial speciation, especially the DNA acquisition capability 
of the bacteria. 

Based on our hypothesis of genetic boundaries and the 
experimental evidence [31], we propose that biological 
species should be non-continuous and non-overlapping 
units of organisms with sufficient similarities among 
members of the same species and radically distinguish- 
ing features from organisms of other species. Bacteria, 
though as essentially non-sexual organisms, should also 
be genetically isolated from one another as species. To 
date, thousands of bacterial strains have been sequenced, 
with a full picture of the genomic divergence among them 
awaiting systematic analyses. We believe that genomic di- 
vergence among bacteria is not continuous and clear-cut 
boundaries ought to exist to demarcate bacteria into 
discrete species. In our Adopt- Adapt Model of bacterial 
speciation [17,32,33], we dissect the process into two over- 
all stages: Adopt, in which bacteria "adopt" exogenous 



DNA from the environment to gain beneficial genetic 
traits, and Adapt, in which bacteria make adjustments to 
"adapt" to the genetic and biological changes brought 
about by the novel genetic material, including genomic 
rearrangements (if the adopted DNA segment was large 
such as larger than 100 kb [34]) and genetic separation 
from the ancestral strain as reflected by the formation 
of the genetic boundaries [31]. According to this model, 
exogenous DNA acquisition is the first step; the present 
study demonstrates that the stochastic lateral DNA ac- 
quisition events may not be destined to take place even 
though the "right" recipient might meet the "right" donor 
at the "right" time (e.g., when the recipient is in bad need 
of certain trait that could be donated by a potential donor 
nearby in the environment). Rather, this event might be 
under strict control, although the molecular mechanisms 
are largely unknown. In this study, we assessed corre- 
lations between genomic sequence divergence and 
transduction or transformation frequency among closely 
related bacteria, anticipating either (i) a linear model with 



Table 3 Relative transduction frequency of Salmonella strains for the five genes 



Strain 


flgL 


treA 


ompC 


hisA 


purG 


mean ± sd 


S. typhimurium LT2 


1 


1 


1 


1 


1 


1 


S. typhimurium 14028 s 


1.17X 10~ 1 


0.60 X10" 1 


0.92 X10" 1 


0.79 x 10" 1 


0.73 x 10" 1 


(8.42 ±2.1 5) x 10" 2 


S. paratyphi B SPB7 


0.42 x 10" 2 


0.24 x 10" 2 


1.03 x10~ 2 


1.36X10" 2 


1.68X10" 3 


(9.46 ± 6.1 1)x 10" 3 


S. agona SARB1 


0.66 X10" 3 


0.11 x 10" 3 


0.61 X10" 3 


1.45 x 10" 3 


0.15 x 10" 3 


(5.94 ± 5.42) x 10" 4 


S. dublin 1TT3627 


4.30 x 10" 3 


1.92 x 10~ 3 


3.64 x10" 3 


3.50 X10" 3 


3.29 x 10" 3 


(3.33 ± 0.87) X10" 3 


S. enteritidis LK5 


2.20 X10" 5 


0.88 X10" 5 


0.00 X10" 5 


1.80 x 10" 5 


1.05 x 10" 5 


(1.19±0.86)x10 -5 


S. paratyphi C RKS4594 


0.93 X10" 3 


0.76 X10" 3 


1.22 X10" 3 


0.13 X10" 3 


0.48 x 10" 3 


(7.02±4.19)x10 -4 


S. paratyphi A ATCC91 50 


1.54X 10" 4 


2.37 x 10" 4 


1.04X10" 4 


0.63 X10" 4 


1.15 x 10" 4 


(1.35 ± 0.66) X10" 4 


S. typhi Ty2 


<1.00x 10~ 6 


<1.00x 10" 6 


<1.00x10" 6 


<1.00x 10" 6 


<1.00x10" 5 


<1.00x 10" 6 


S. arizonae SARC5 


8.15 x 10" 4 


3.08 x10~ 4 


5.53 X10" 4 


1.17 x 10" 4 


3.46 x 10" 4 


(4.28 ± 2.66) x 10" 4 


S. bongori S ARC1 2 


<1.00x 10~ 6 


<1.00x 10~ 6 


<1.00x10" 6 


<1.00x 10" 6 


<1.00x10" 6 


<1.00x 10" 6 



Note: The relative transduction frequency listed in the table refers to relative transduction rate compared to LT2. 
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Figure 2 Transduction frequencies against levels of relatedness 
among the Salmonella strains. 



falling transduction or transformation frequency in cor- 
relation with increasing genomic divergence; or (ii) a 
non-linear model; we observed the latter. 

To interpret these observations, we may divide the ex- 
ogenous DNA acquisition process into two steps. In the 
first step, exogenous DNA enters the cell and escapes 
degradation by restriction enzymes or other defense sys- 
tems. The next step includes hetero duplex molecule 
formation, replication and segregation to daughter cells. 
This second step is controlled negatively by mismatch 
repair (MMR) proteins, such as MutS and MutL, and 
positively by the induction of SOS system. As MMR and 
SOS genes are highly conserved in bacteria, they are 
supposed to function equally in all tested Salmonella 
strains. Consequently, the second step here should mainly 
be determined by the sequence divergence between the 
donor and recipient. It is therefore reasonable to presume 
that the non-linear correlation between recombination 
efficiency and genomic divergence should be due to 
differences in the first step, i.e., entry of exogenous 
DNA into the bacterial cell. As transduction and trans- 
formation use different mechanisms to bring DNA into 
the cell across the cell wall and membrane, the plasmid 
transformation may not effectively simulate the first step of 
P22-mediated homologous recombination. Nevertheless, 
when we divided the transduction frequency by the trans- 
formation frequency in a given strain, the "normalized" re- 
sults fitted the linear model slightly better (calculation not 
shown). Taken together, our results demonstrate that the 
linear model may hold true only when recombination effi- 
ciency is determined exclusively by sequence divergence, 
although there seems no such a case in the real world. As 
this study focused on only DNA transfer events, systematic 
elucidation of the bacterial speciation mechanisms requires 
investigation of also other strain differences such as level 
of gene expression for competency/DNA recombination, 



which is under study in our laboratory. Since transform- 
ation and transduction are the two main mechanisms for 
DNA uptake by bacteria and even very closely related bac- 
teria may exhibit distinct performance in recombination 
as demonstrated in this study, it is possible that subpopu- 
lations or even individual cells of a bacterial species might 
have begun to form genetic barriers against exogenous 
DNA to stabilize the bio-system within the cell or a 
subpopulation or to continue the speciation process to 
become a nascent species. 

Conclusion 

DNA acquisition efficiency is variable among closely re- 
lated bacteria, suggesting that bacterial speciation involves 
highly complex processes affected not only by whether 
beneficial exogenous DNA may exist in the environment 
but also the "readiness" of the bacteria to accept it. 

Methods 

Bacterial and phage strains 

The Salmonella strains used in this study are listed in 
Table 1; more detailed information on these bacterial 
strains can be found at the Salmonella Genetic Stock 
Center (SGSC; http://www.ucalgary.ca/~kesander/). The 
bacterial strains were cultured at 37°C and the bac- 
teriophage P22 strain (HT105/1 int-201) used in the 
transduction experiments was routinely grown on S. 
typhimurium LT2. 

Phage P22-mediated transduction 

P22-mediated transduction was used to mediate homolo- 
gous recombination between donor and recipient strains 
as previously described [35,36] . The primary donor strains 
were derivatives of S. typhimurium LT2 [37,38]. Five such 
derivative strains were used, having TnlO (with tetracyc- 
line resistance) inserted within gene flgL, treA, ompQ hisA 
or purG. After P22 infection of the secondary donor 
strains (S. gallinarum or S. pullorum strains used as 
donors of the Tn70-inserted genes to the S. gallinarum 
or S. pullorum recipient strains; see Results), the phage 
lysate was used to infect recipient strains. Construction 
of the secondary donor strains from the TniO-inserted 
LT2 derivatives and preparation of phage lysate were 
done as described [36]. 

As P22 can only infect the bacteria that express the 
012 antigen, we first transformed the 012"recipient 
strains by a cosmid, pPR1347, that carries the 012 anti- 
gen gene [39]; details of the experimental procedure 
can be found in literature [36]. Briefly, the E. coli strain 
that carries pPR1347 was cultured in LB broth over- 
night. The cosmid was extracted by alkaline lysis method 
and was transferred into S. paratyphi C RKS5478, S. arizo- 
nae SARC5 and S. bongori SARC12; all of the three 



Bao et al. BMC Microbiology 2014, 14:1 57 
http://www.biomedcentral.com/1471 -21 80/1 4/1 57 



Page 7 of 9 



£2 o 



U 3 



Groups 




12 12 
Groups Groups 

Figure 3 Comparisons of recombination efficiency reflected by transduction frequencies, with the points representing the numbers of 
transductants. A, Transduction of ompD159 between S. gallinarum and S. pullorum (see text for more details). Groups: 1, S. pullorum recipients and S. 
pullorum DNA; 2, S. pullorum recipients and S. gallinarum DNA; 3, S. gallinarum recipients and S. gallinarum DNA; 4, S. gallinarum recipients and S. 
pullorum DNA. Note that group 3 has the largest numbers of recipients (groups 2 vs 3, p = 0.007; groups 3 vs 4, p = 0.000), reflecting significantly higher 
recombination efficiency of S. gallinarum DNA into the 5. gallinarum genome (group 3) than 5. gallinarum DNA to the 5. pullorum genome (group 2) or 
than S. pullorum DNA to the 5. gallinarum genome (group 4). Group 1 had a similar tendency as group 3, although the differences were not statistically 
significant (groups 1 vs 2, p = 0.079; groups 1 vs 4, p = 0.071). B, Transduction with S. pullorumas the recipient and S. pullorum or 5. gallinarumas the 
donor. Groups: 1, S. pullorum recipients with DNA fromS. gallinarum 287/91; 2, S. pullorum recipients with DNA from 5. pullorum RKS5078 (groups 1 vs 2, 
p = 0.009). C, Transduction with S. pullorum or 5. gallinarumstralns as the recipient and donor DNA from 5. gallinarum. Groups: 1, S. pullorum recipients 
with DNA from S. gallinarum strain 287/91; 2, S. gallinarum recipients with DNA from 5. gallinarum strain 287/91 (groups 1 vs 2, p = 0.000). 



Table 4 Frequencies of transformation among the 
recipient strains for the same donor DNA 



Strain Transformation frequency 



S. typhimurium LT2 


(4.30 ± 1.15) X10" 4 


5. typhimurium 14028S 


(5.08 ± 0.72) X10" 6 


S. paratyphi B SPB7 


(4.1 0 ± 1.40) x 10" 4 


S. agona SARB1 


(2.25 ± 0.90) X10" 6 


S. cfub//n m3627 


(2.46 ± 0.91) X10" 4 


S. enteritids LK5 


(8.33 ± 1.44) x 10" 7 


S. paratyphi C RKS4594 


(1.44±0.01)x 10" 4 


S.paratyphi A ATCC9150 


(2.93 ± 1.05) X10" 4 


S. typhi Jy2 


(8.15 ± 1.43) X10" 4 


S. arizonae SARC5 


(1.32±0.51)x10~ 3 


S. bongo/7 SARC 12 


(1.58 ± 0.45) x 10" 3 



bacterial strains expressed 012 after the transformation 
by pPR1347. 

The transduction frequency was determined as fol- 
lows. The recipient strains were cultured in LB broth till 
OD600 reaching 2.0. Then the bacteria were mixed with 
P22 phage lysate at an MOI (multiplicity of infection) of 
0.01, plated onto LB medium containing tetracycline 
(15 (ig/ml) and incubated at 37°C for 18 h. Because all re- 
cipient strains were sensitive to tetracycline, only clones 
that have gone through successful homologous recombin- 
ation would be able to grow on the tetracycline-containing 
plate. The relative transduction frequency of each recipient 
strain was calculated by dividing the number of its trans- 
ductants by the number of LT2 transductants transduced 
by the same lot and amount of the phage lysate. Each 
experiment was repeated five times. 
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Transformation 

The plasmid pQE30, which carries ampicillin resistance, 
was used for transformation. The recipient bacteria were 
cultured in LB broth until OD600 reached 0.4. Compe- 
tence was induced by thawing the culture on ice for 
1 min. A volume of 80 \A competent cells were mixed 
with 2 [A plasmid pQE30 and electroporated at 25 (iF 
capacitance, 2.5 kv and 200 ohm resistance, followed 
immediately by the addition of 320 (il fresh LB medium. 
The bacteria were incubated on 37°C with agitation of 
85 rpm for 1 hour. Then the bacteria were spread onto 
LB agar plates containing ampicillin (100 (ig/ml) and in- 
cubated at 37°C overnight. As all recipient strains used 
were sensitive to ampicillin, only successful transformants 
could grow on ampicillin plates. The transformation 
frequency was calculated by dividing the number of 
transformants by the number of competent cells of the 
same volume grown on pure LB plate. Each experiment 
was repeated five times. 

Calculation of sequence divergence 

The genome sequences for comparison were downloaded 
from Genbank database, with the accession numbers 
being listed in Table 1. The genome sequence of S. 
agona SARB1, S. dublin TYT3627, S. enteritidis LK5 are 
not available in public database yet, so we use the pub- 
lished strains of the same serovar instead (see Table 1). 
Because S. bongori SARC12 doesn't have any published 
relative of the same serovar, we sequenced its genome 
by using SOLiD™ 3 0 sequencer. Briefly, we sheared the 
bacterial genomic DNA into fragments of 2 ~ 4 kb in 
size and sequenced both ends of the fragments follow- 
ing SOLiD™ 3 0 2 x 50 bp mate-pair sequencing protocol. 
The raw reads were assembled into draft genome by using 
the software SOLiD™ System de novo Accessory Tools 2.0 
(http://solidsoftwaretools.com/gf/project/denovo/). The gaps 
between the obtained contigs were closed by PCR. 

We calculated the sequence divergence between the 
compared strains by two approaches. The first approach 
was based on five genes only. The nucleotide sequences 
of flgL, treA, ompC, hisA and purG were extracted from 
the complete genomes, aligned by CLUSTALW, and 
then concatenated together. The alignment was input 
into MEGA5 for constructing the phylogenetic tree. 
The Jukes-Cantor substitution model was adopted and 
the neighbor-joining algorithm was implemented. The 
second approach was based on all conserved ortholo- 
gous genes. We chose LT2 as the reference genome and 
searched its protein-coding sequences against other 
strains by BLAST. The threshold of the BLAST search 
was as follows: the e-value must be smaller than 1 x 10" 10 , 
the identity of the amino acid must be larger than 80% 
and the alignment length must be longer than 80% of the 
gene length. The conserved genes among all strains were 



then aligned by CLUSTALW and were concatenated to- 
gether. Finally the alignment was input into MEGA5 and 
the construction of the phylogenetic tree was built with 
the same parameter to the first approach. 

Availability of supporting data 

The genome sequence of S. bongori SARC12 is available 
the Genbank database under the accession number 
CP006692. The phylogenetic trees supporting the re- 
sults of this article are available in the TreeBASE re- 
pository, http://purl.org/phylo/treebase/phylows/study/ 
TB2:S15939?x-access-code=26f7bb0963bf24cff2e944fe3 
065c5ec&format=html. 
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